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Abstract 

It is shown that the calculation of the magnetic moment of a 
Friedel-Anderson impurity in mean-field theory is unreliable. A class 
of approximate solutions, which contains the mean-field solution as 
an element, is expressed in rotated Hilbert space and optimized. The 
optimal state has considerably lower energy than the mean field so- 
lution and requires almost twice the Coulomb exchange U to become 
magnetic. Since most moment calculations of magnetic impurities, 
for example the spin-density-functional theory, use the mean-field ap- 
proximation the resulting magnetic moments have to be critically re- 
examened. 

PACS: 75.20.Hr, 

The properties of magnetic impurities in a metal is one of the most inten- 
sively studied problems in solid state physics. Although some of the exper- 
imental anomalies were already discovered in the 1930's, it is still a subject 
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of great interest. The work of Friedel PQ and Anderson {2| laid the founda- 
tion to understand why some transition metal impurities form a magnetic 
local moment, while others don't. They considered a host with an s-band in 
which a transition metal atom is dissolved. The s-electrons can hop onto the 
d- impurity via the hopping matrix element V s d- The ten- fold degeneracy of a 
real d-impurity is simplified and reduced to a two-fold degeneracy for spin up 
and spin down. If both states are occupied they repel each other due to the 
Coulomb exchange energy. This yields the Friedel- Anderson Hamiltonian 

N N 

Hp A = ^{^^C^C^ + ^dC^ + ^Kdl^Kc^ + C^^I + ^d+^d- (!) 
<j v=l u=l 

Here a finite s-band with N states is used. The c* CT and the d* are the 
creation operators of the (free) s-electrons and the d-impurity. The (i*-states 
are assumed to be orthogonal to the s-states c* (in the following I denote 
single electron states by their creation operator). 

In the limit of V 8 d = and E d < ep,E d + U > e F the d-impurity is 
magnetic. Anderson concluded that the magnetic moment survives for small 
but finite V s d and derived the criteria for a magnetic state and the size of 
the moment in a mean-field approximation. He found a magnetic state if 
the product of Ug^ > 1 where g d is the additional density of states of the 
d-resonance. 

Kondo |2] brought a new twist into the magnetic impurity problem when 
he showed that multiple scattering of conduction electrons by a magnetic im- 
purity yields a divergent contribution to the resistance in perturbation the- 
ory. In the following three decades a large number of sophisticated methods 
were applied to better understand and solve the Kondo and Friedel- Anderson 
model, and it was shown that at zero temperature a Friedel- Anderson impu- 
rity is in a non-magnetic singlet state. However, above the Kondo tempera- 
ture the impurity shows a magnetic moment, and there is a great interest in 
the size of this moment. 

There is a large body of research in which the magnetic moment of im- 
purities is calculated jlj, [5], jSj, [Z|, [H] Generally spin-density-functional 
theory is used for this task. Within this theory the electronic structure of 
the host and the impurity is calculated from first principles without any 
adjustable parameters. In particular the strength of the Coulomb and ex- 
change interaction are obtained from first principles. However, in the final 
step the mean-field method is applied to obtain the local magnetic moment. 
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Although this is a zero-temperature calculation (where the impurity should 
be in the Kondo singlet state) it is generally argued that such a calculation 
yields the magnetic moment above the Kondo temperature (which, at lower 
temperatures, is hidden in the singlet state). 

In this letter I will show that the mean-field result for the magnetic mo- 
ment of impurities is not reliable. By rewriting the mean-field solution in a 
rotated basis and optimizing the solution I obtain solutions which are much 
lower in energy, require a much larger critical U for the formation of a mo- 
ment and yield smaller moments. And this despite the fact that the improved 
solution has the same structure (in the rotated basis) as the mean-field solu- 
tion. Since there is a large body of spin-density-functional theory calculations 
for magnetic impurities, a reevaluation of this method might be required. 

I start with Anderson's (potentially) magnetic state which he obtained 
as a mean-field solution. Anderson replaced the Hamiltonian Hfa by 

H mf = H F+ + H F _ - U(n d+ )(n d _) (2) 

N 

H Fa = E u e u cl a c ua + ^2 Vsd{v) [d* a c ua + c* CT rf CT ] + E d>a d* a d a (3) 

u=l 

where (n d+ ) and (rid-) are the average occupation numbers of the states d* + 
and d*_ and E d:(7 = (E d + U (n d - a )). The solution of the mean-field method 
requires the diagonalization of two Friedel resonance Hamiltonians H Fa with 
self consistent values for (n d+ ) and (n d J). This straight-forward numerical 
calculation yields the mean-field ground-state energy E m f and the magnetic 
moment fj, m f. The energy of the bare magnetic state E b rn _ is subtracted from 
E m f where 

n 

E b . m .=2j2^ + E d -e n (4) 

is the ground-state energy for V sd = and E d < e F , E d + U > e F . 

For the numerical calculation an s-band with a constant density of states 
is used, ranging from —1 to +1. This band is divided into iV = 48 equal 
cell. Each s-sub-band is half filled, i.e. the number of occupied states in 
each spin sub-band is n = AT/2. In Fig.l the numerical results for E m t — E b , m , 
are plotted for |V^| = 0.05. The Coulomb repulsion U is varied between 
0.2 and 1.2. Together with the Coulomb repulsion the <i*-state energy E d is 
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varied so that and (E^ + U) lie symmetrically about the Fermi energy, 
i.e.E d = -\U. 

In the mean-field calculation the impurity is non-magnetic for U < U cr ~ 
0.275. For U > U cr the spin up and down sub-bands split. The resulting 
magnetic moments are plotted in Fig. 2 (curve with circles). 
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Fig.l: A comparison between the ground-state energies of the mean-field 
calculation and the AFR method (see text below). 



Since the mean-field solution is the product of two n-electron states of 
the two Friedel Hamiltonians Fp + and Hp_ we consider these solutions of 
the Friedel Hamiltonian (JHJ) in some detail. As shown in ref. ^0] the 
exact ground state of Hp with n (spinless) electrons can be written in the 
form 

n-l 

^ F =[A'a* + B'd*}Y[a*% (5) 

i=l 
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Here <3>o is the vacuum state and (Iq is cl sister state to d* which is built 
from the states of the s-band 

JV 
v=\ 

Ref. ^Uj shows how to calculate the coefficients oP v from the parameters of 
the Hamiltonian Hp and the occupation number n. The a* are orthogonal to 
cIq and to each other and their (N — 1) sub-matrix of the s-band Hamiltonian 
H = ^2e u n u is diagonal (see equ. |SJ). The states a* are uniquely determined 
from the state Oq. Their form is 

N 

a* = <*i< ( 7 ) 

v=\ 

The a* (1 < i < N — 1) together with represent a new basis. 

In this new basis the free electron Hamiltonian H Q = J2u=i £ u c t°v takes 
the form 

N-l N-l 

H = Y, E (*) a >i + E (0) ao«o + Yl V fr (*) l a o a i + a >o] ( 8 ) 
1=1 1=1 

In the Hamiltonian (jSJ) the Og-state represents an artificial resonance state. 
I will call it in honor of Friedel an Artificial Friedel Resonance state (AFR 
state). It is a sister state to the state d*. 

The full (spin independent) Friedel Hamiltonian can be written as 

N-l 

H F = ^E (i) a*ai + E (0) a* a + E d d*d + V s a d {0) [d*a + a*d] (9) 

i=l 

N-l N-l 

i=l i=l 

where 
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In the Hamiltonian Q the d*-state and the localized a^-state are on 
equal footing. The second line in equ. Q yields the hopping between a* 
and d* (first term) and a* and (second term). For the state (A'clq + B'd*) 
the individual hopping matrix elements cancel each other, making typ r the 
ground state. 

In the next step the mean-field solution is rewritten in the AFR-form 
of the Friedel ground state. Since the Hamiltonian consists of a Friedel 
Hamiltonian for each spin the mean-field state is the product of two states 
of the form of equ.(|H|l. Therefore this mean-field state (the exact solution of 
the mean-field Hamiltonian) can be written as 

n-l 

= [A-a*_ + B„d*_ ] [A+a* 0+ + B + d* + ] J] a* a % 

n-l 

= [Aa* _a* 0+ + Bd*_a* 0+ + Ca* ^d* + + Dd*_ d* + ] JJ a* ia % (11) 

a,i=l 

= A$ a + B^ B + C^c + D^ D 

where 

A 2 + + Bl = 1 , A 2 _ + B 2 _ = 1 

A = A + A^ , B = A + B_ (12) 
C = A^B + , D = B + B_ 

Each of the four states ^ a, ^b, and ^ d is normalized, and they are all 
orthogonal to each other. In the magnetic solution one has A + ^ A^ and 
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B + 7^ B_. Also the two rotated bases {a^ + , a* + } and {a^— ? a i-} are different 
in the magnetic state. 

So far the many electron state in equ. (jSJ) is the mean-field solution. This 
state consists of an electron background n™7=i a ia^o multiplied with the sum 
of four two-electron states, consisting of the combinations [aQ_a,Q + , d^a,Q + , al_d* + , d*_d\\ 
which have S z = 0. The mean-field wave function opens a new playing field 
for variation to find the optimal state. One can optimize the coefficients A, 
B, C and D while dropping the individual normalization conditions (jl2j) and 
replacing them by 

A 2 + B 2 + C 2 + D 2 = 1 (13) 

Far more important one can optimize the states a$ + and ao_. 

For this purpose the Hamiltonian Hp a is expressed in the bases {a^ + , a* + } 
and {aQ„,<2*_} . One obtains for the expectation value of the ground-state 
energy E 

E = A 2 [£L (0) + E + (0)] + B 2 [£L (0) + E d \ + C 2 [E+ (0) + E d \ + D 2 [2E d + U\ 

n-l 

+ 2(AB + CD) V- d (0) + 2 (AC + BD) V+ (0) + £ E a (i) (14) 

For a given set of states {aQ ± ,a* ± } the energy £"0 in eq. (jlljl depends 
on the coefficients A, B, C and D. One obtains the lowest energy by vary- 
ing Eq with respect to these coefficients. This yields a 4x4 matrix for the 
coefficients vector (A, B, C, D). The lowest eigenvalue gives the energy ex- 
pectation value, and its eigenvector gives the coefficients. The resulting state 
I denote as the magnetic state ^afr and the solution as the AFR solution. 

The central part of the numerical calculation is the variation of the states 
<2q + and <Zq_ until the absolute minium of the energy is reached. 

As in the mean-field theory the numerical calculation itself determines 
whether the lowest state possesses a magnetic moment or not. If the solution 
is magnetic then and approach different states and the coefficients 
B and C have different values. The resulting magnetic moment is defined as 
the difference in the occupation of the d*^ and d*_ states, i.e. /1 = B 2 — C 2 . 

In Fig.l the energy expectation value E of the optimal magnetic state 
^ afr is plotted as the curve with the triangles (again the same energy E b m ^ 
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has been subtracted). The new ground-state energy lies considerably below 
the mean-field energy. 

In Fig. 2 the resulting magnetic moments that one obtains with the mean- 
field approximation and with the new method are plotted. One recognizes 
that the new solution suppresses the magnetic moment up to a considerably 
larger value of U cr 0.46. This is almost twice the value of the mean-field 
theory. 




Fig. 2: The magnetic moment as a function of the Coulomb energy U, using 
the mean-field solution and the AFR-method of the present calculation. 



How do we have to interpret the fact that the AFR solution suppresses the 
magnetic moment up to a much larger critical Coulomb exchange interaction? 
Since this state has a lower energy expectation value, does this mean that 
its magnetic moment is more reliable? The author prefers a somewhat dif- 
ferent interpretation. The Friedel-Anderson impurity does not like a broken 
symmetry. The mean-field approach does not give the multi-electron state 
any wiggling room. Only the values for E^ a can be altered with increasing 
U. The AFR solution on the other hand possesses a lot more flexibility since 
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the AFR states can adjust. Therefore the symmetric multi-electron state 
survives to a considerably larger Coulomb exchange interaction. It might be 
that neither state yields the right magnetic moment for the impurity. The 
present calculation raises serious questions about the mean-field approach. 
This might also apply to the spin- density-functional theory (SDFT) for mag- 
netic impurities. This theory is a very complex theory and it is difficult to 
judge from the outside all the intricacies. It should yield the correct charge 
and spin densities for the correct functional. But in the final step the ma- 
jority of SDFT calculations use a two-spin-fluid model where each electron 
fluid adjusts in the (mean) field of the other. 

To conclude, in this paper an approximate solution for the Friedel- Anderson 
impurity is expressed in a rotated Hilbert space {clq ± , a i±}- Its center piece 
are two artificial resonance states ag+' a o- f° r the spin up and down s- 
electrons. They determine uniquely the remaining bases {a* + },{a*_}. The 
AFR states are combined with the d-electrons for spin up and down d* +1 d*_ 
into two-electron states of total S z = 0, i.e. [Aa>Q_aQ + + Bd*_a^ + + Ca^_d* + + Dd*_d* + ] . 
Then the {n — 1) lowest states of the two (N — 1) bases {a*±} are occupied 

71-1 

yielding the s-electron background Y\ a ia^o- The compositions of the AFR 

i=l,a 

states <2q + , <2q_ are calculated by numerical variation which rotates the s- 
electron bases in Hilbert space. 

The energy of the resulting state lies clearly below the mean-field solution. 
The critical value of the Coulomb exchange energy U cr for the formation of a 
magnetic moment is almost twice as large as in the mean-field solution. Since 
in many calculations of the magnetic moment of impurities the mean-field 
approximation is used one has to reevaluate the resulting moments. This may 
also apply to the impurity calculations which use the spin-density-functional 
theory because in the majority of these calculations the mean-field theory is 
used in the final analysis. 

Since the ground state of the Friedel- Anderson impurity is a singlet state 
one might suspect that the structure of the new solution with the lower 
energy and smaller magnetic moment is somewhat closer to the singlet state 
than the mean-field solution. This is not the case. Both the mean-field and 
the present solution are in a symmetric state for small U ; both show a similar 
asymmetry between spin up and down in the magnetic state. The mean-field 
solution belongs to the same class of wave functions as the here presented 
one (which are given by the general form of equ. (jlljO . 
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